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ABSTRACT 

This  paper  represents  aui  attempt  to  obtain  useable  approximations, 
for  two  special  cases,  to  the  scattering  matrix  for  an  inhomogeneous,  aniso- 
tropic medium  of  finite  thickness  "bounded  ahove  and  helow  hy  semi-infinite 
homogeneous,  isotropic  media.  The  ultimate  purpose  of  this  attack  will  "be 
to  obtain  the  transmission  and.   reflection  coefficients  for  a  plane  wave  in- 
cident at  an  arbitrary  angle  on  an  ionospheric  layer  with  constant  e  and  yi, 
in  a  Tiniform  earth's  magnetic  field  as  a  functional  of  both  the  charge  dis- 
tribution function  N(z)  and  the  collision  frequency   »>  (z). 

First,  the  scattering  matrix  for  an  anisotropic  multilayer  is  ob- 
tained. This  will  be  a  good  approximation  to  the  continuously  varying  medium 
at  low  frequencies  where  each  layer  is  but  a  fraction  of  a  wavelength  in  thick- 
ness. 

In  the  high  frequency  case,  where  a  solution  in  series  form  has  al- 
ready been  obtained  by  Bremmer,  it  is  found  that  a  considerable  simplification 
can  be  effected  by  expanding  the  conductivity  matrix  in  inverse  powers  of 
frequency. 
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1,   INTRODUCTION 

The  use  of  electromagnetic  waves  for  determining  the  structure  and  pro- 
perties of  the  ionosphere  has  stimulated  interest  in  investigating  the  behavior 
of  these  waves  in  a  horizontally  stratified  ionized  medium  in  the  presence  of  a 
fixed  external  magnetic  field.  Such  a  medium  is  a  reasonahle  approximation  to  the 
ionosphere.  The  problem  is  most  readily  treated  in  the  range  of  high  frequencies, 
where  recourse  can  he  had  to  verious  modifications  and  improvements  of  the  W.Z.E. 
approximations  or,  alternatively,  to  the  series  expansion  technique.  The  lowest 
order  of  approximation  in  the  latter  case  is  the  geometrical  optics  approximation  [l]. 

When  the  frequency  of  the  incident  electromagnetic  wave  is  extremely  low 
(below  100  E.G.),  neither  the  series  expansion  technique  nor  the  W.K.B.  approxima- 
tions are  applicable.   In  this  frequency  range,  the  anisotropic  multilayer  solution 
of  Section  2  may  prove  useful. 

This  method  consists  of  approximating  the  ionosphere  by  a  piecewise  con- 
stant medium  composed  of  a  number  of  homogeneous,  anisotropic  slabs,  each  of  which 
is  but  a  fraction  of  a  wavelength  in  thickness.  The  cese  of  an  isotropic  multi- 
layer hfls  been  treated  by  many  writers  \Zj   and  the  method  used  here  is  basically  a 
generalization  of  the  procedure  used  in  that  case. 

For  a  given  Incident  plane  wave,  the  electromagnetic  field  in  each  Blab 
is  known  to  consist  of  a  linear  combination  of  four  elementary  plane  waves,  two 
upgolng  and  two  downcoming.  The  coefficients  of  these  elementary  plane  waves  are 
related  by  the  contiiiuity  conditions  on  the  electric  and  magnetic  field  across  every 
Interface  separating  a  pair  of  slabs.   In  this  manner,  it  is  possible,  by  ordinary 
algebraic  manipulation  to  determine  the  amplitudes  of  the  outgoing  plane  waves  in 
terms  of  the  ainplitudes  of  the  incoming  plane  waves.  This  connection  is  given  by 
a  matrix  known  as  the  ^'scattering  matrix". 

As  the  number  of  slabs  in  the  multilayer  is  made  infinite,  the  thickness 
of  each  slab  going  to  zero,  the  multilayer  becomes,  in  the  limit,  a  continuous 
medium.  The  multilayer  solution  can  be  snown  \.o   re&uce,  in  this  limit,  to  an  iterate* 
solution  obtained  by  Bremmer  [3]  for  the  electromagnetic  field  in  a  continuously 

[1]   See,  for  example,  H.G.  Booker,  Phil.  Trans.  Hoy.  Soc.  237,  ^11,  1938 

[2j   See,  for  example,  E.K.Luneberg,  The  Propagrtion  of  Electromagnetic  Plane  Waves 
in  Plane  Parallel  layers.  Besearch  Report  No. 172-3, W.S.C.  Ms.th.Res.Grp.  19^7 

[3]   H. Bremmer,  Terrestrial  Radio  Wav«s,  pp. 301-336  Elsevier  Publishing  Co.,  N.Y., 
1950. 
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varying  anisotropic  medium.  While  the  latter  solution  is  exact,  it  is  in  series 
form,  and  the  series  is  rapidly  convergent  only  at  high  frequencies.  For  practi- 
cal purposes  it  can  therefore  he  used  only  in  this  range.  The  multilayer  solution, 
on  the  other  hand,  is  applicable  only  in  the  low  frequency  range,  since  each  layer 
of  the  mulTuilayer  must  he  considerably  smaller  than  a  wavelength  in  thickness  and, 
for  short  wavelengths,  the  number  of  layers  would  have  to  be  too  large  for  practical 
computation* 

In  Section  3  a  simple  formula  is  derived  for  the  first  few  terms  in  a 
series  expression  for  the  scattering  matrix  as  explicit  functions  of  the  ionospheric 
parptneters  for  the  case  in  which  the  frequency  of  the  incident  wave  is  greater  than 
both  the  maximum  critical  frequency  and  the  gyromagnetic  frequency  of  the  iono- 
spheric layer.  The  starting  point  of  this  discussion  is  the  Bremmer  solution,  which, 
although  rapidly  convergent,  is  given  in  terms  of  complicated  functions  of  the 
solutions  of  a  quartic  equation.  By  expanding  the  quantities  appearing  in  the  wave 
eqiiation  for  the  electric  field  in  inverse  powers  of  the  frequency  of  the  incident 
wave,  a  technique  used  by  Ginsburg  f^J  ,  the  algebraic  complexity  of  the  Bremmer 
solution  is  greatly  simplified. 

The  simultaneous  use  of  the  series  expansion  technique  and  the  iterated 
solution  technique  has  a  further  substantial  advantage.  The  formula  obtained  in- 
cludes the  effect  of  the  collision  frequency  on  the  electromagnetic  field  explicitly, 
whereas  this  is  ordinarily  too  complicated  to  be  handled  by  either  the  iterated 
solution  technique  or  the  series  expension  technique  alone. 

The  technique  used  in  section  3  to  treat  the  high  frequency  case  may 
also,  with  suitable  modifications,  be  used  for  the  case  of  low  charge  density. 
This  problem  will  form  the  subject  of  another  report.  Only  with  this  approxi- 
mation, along  with  the  two  approximations  treated  in  this  report  can  nunerical 
computations  on  actual  ionospheric  models  effectively  be  performed.  For  this 
reason  no  numerical  work  is  included  in  this  papeiv 

a 


) 


LU]  V,  L.  Glnsburg,  Journal  of  Physics,  Vol.  VII  No.  6,  289  (19^5) 
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2.      THE  GEN5BAL  A^T ISOTROPIC  MULTlUrEB 

In  this  section  we  will  consider  an  obliquely  incident  plane  wave 
in  a  strfitifled  anisotropic  medium  which  consists  of  N  +  2  homogeneous 
strata.     The  zeroth  stratum,  which  extends  from  z=-ootoz=0  and  the 

St 

N  +  1   stratum,  which  extends  from  z  =  L  to  z  =  +oo  will  both  be  assumed 
to  be  isotropic;  the  other  N  strata  will  be  restricted  only  in  that  the 
permeability  p.   (u=l,..,,N)  will  be  requii"ed  to  be  scalar.  Hence  e 
and  C        may  be  taken  to  be  arbitrary  tensor  quantities,  although  the 
most  useful  case  will  be  that  in  which  p.    and  e    are  both  scalars 
and  the  conductivity  tensor  0'^"'  is  that  resulting  from  a  free  charge  dis- 
tribution in  a  uniform,  earth's-  magnetic  field.   If  N  is  then  taken  sufficient- 
ly large,  this  special  case  will  yield  a  good  approximation  to  t^he  physical 
sit-uation  of  a  plane  wave  obliquely  incident  on  the  ionosphere. 

For  the  sake  of  simplicity,  the  x-  and  y-  coordinate  axes  will  be  so 
chosen  that  the  propagation  vector  in  the  zeroth  medium  (and  hence  in  all 
media)  lies  in  the  x-z  plane,  the  z-axis  having  alreaay  been  chosen  per- 
pendicular to  the  strata. 


^  =  ^N-l 

^=  «^1 
Z  =  0 


/  X  N  ^-  1   medium   X  / 


th  _,. 
V       medium 


//  /  / 

Fig.  1 


0'"  medium  y  /  / 


Under  these  conditions,  the  electric  field  E    in  the  v       medium 
can  be  expressed  as  a  sum  of  the  form* 

(1)   E(''\x,z.t)=i:  C^^^-'V  ^"'^^^^n!^ 

/  »         n=l  _/  \ 

where  C    is  a  complex  amplitude,  A    is  a  complex  polarization  vector 
n  n 

which  depends  only  on  the  properties  of  the  medium  as  does  the  z-component-^f 


•  The  justif Icetion  for  assuming  an  expansion  of  the  form  (1)  for  E  will 
be  given  in  Section  2. A.   It  is  also  disciassed  in  H.  Bremmer,  Terrestrial 
Radio  VJaves,  pp.  31-0-316, 
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the  propagation,  vector  h   .  The  value  of  the  x-component  of  the  propagation 

vector  k  is  known  to  "be  independent  of  the  medium  and  depends  only  iipon  the 

direction  of  incidence  and  the  frequency  of  the  incident  wave. 

Our  first  problem  (section  Ik)   will  he  to  determine  A   and  h    in 

n      n 

each  of  the  media  as  a  function  of  the  physical  properties  of  that  medium. 
Then,  in  section  IB,  by  applying  continuity  conditions  across  the  "boundary 
separating  the  v       medium  from  the  u-1  ,  we  shall  he  ahle  to  determine  the 
amplitudes  C    as  a  linear  comhination  of  the  amplitudes  C     in  the  slat 
below.  3y  repeating  this  process,  we  can  determine  the  amplitudes  in  any 
slab.  Thfen,  using  equation  (l),  we  can  specify  the  electric  field  at  any 
point  in  the  medium.  From  this,  the  values  of  the  components  of  the  magnetic 
field  can  be  computed  from  Maxwell's  equations,  if  desired. 
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^  "  Determination  of  the  Propagation  and  Polarization  Vectors. 

If  we  assTioe     fL     and     ?     to  te  arbitrary  tensors  which  do  not  vary 

with  time  and  if  we  assume  that  I,  H  and  J  have  the  harmonic   time  de- 

— iwt 
pendence  e  ,  we  can  write  Maxwell's  equations  as 

(2)  ^  "  +i«J'lS 
7x1"  =  -ittieff  +  J 

where   pH  and  Ts^are  the  vectors  resulting  from  the  matrix  maltiplication 
of  ^  and  H  and  €"  and  E  respectively.   If,  now,  we  restrict  ?  to  he  a 
scalar  quantity  designated,  simply,  hy   ^,  we  can  eliminate  H  f rom  (2) 
to  get 

(3)  y-^^vxl^  m^^,-^^  i  ^J  ^ 

We  now  replace  J  by  ^,  where  ^  is  an  arbitrary  conductivity  tensor  and 
rewrite  (3)  as 

where  we  have  abbreviated  the  right  hand  side  by  writing 
(4b)       p=  uj  2  ^(  7  +  ^^)  . 

Now,  in  each  of  the  media,  n  is  constant,  so  that  the  left  hand  side 
of  {U&)   becomes  ^  x  V  x  I,      Substituting  for  E  a  plane  wave  solution 

E=Ae^  [- '"t+kx-^hz] 
into  (U),  we  get,  writing  out  the  equations  fully, 

The  set  of  linear  equations  (5),  being  homogeneous,  will  possess  a  non-trivial 
solution  if  and  only  if  the  determinant  of  the  coefficients  vanishes.   Setting 
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this  determinant  equal  to  zero  yields  the  following  quartic  eq\iation  in  h 
for  a  fixed  value  of  k. 

.(">  „(")  pt") 

xy        xz 


(u)  ,  p(.)j 
XX     zz 


xx 

M 

yx 

(a) 

zx 


77 
zy 


7z 

(w) 
zz 


} 


=  0 


.th 


Thus,  there  will,  In  general,  he  four  values  of  h  in  the  v       layer  for 
which  plane  wave  solutions  of  (U)  with  a  given  value  of  k  exist.  We  will 
designate  them  hy  h    (n=l,2,3,^)  and.  the  corresponding  polarization  vectors, 
obtained  hy  substituting  these  values  for  h^  '  into  (5),  by  A   .  The  most 
general  plane  wave  solution  of  the  wave  equation  in  the  u   medium  will  there- 
fore be  a  linear  combination  of  the  four  elementsury  solutions  we  have  just 
obtained  multiplied  by  four  amplitude  coefficients  C^"  ,  thereby  justifying 


the  form  for  E  assumed  in  equation  (l). 

As  stated  before,  the  polarization  vectors  A    can  be  determined  (to 
'  n 

within  an  arbitrary  normalization  factor)  by  solving  equations  (5);  that  is, 
the  ratio  of  any  two  components  of  A    can  be  determined  from  (5).  Thus, 


we  have 


Xv) 
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where  m    la  an  artitrary  normalization  factor. 


n 


.(") 


It   1b   Interesting  to  note   that  the  quartic  equation   (6)  for  h         reducea 
to  a  'bi-quadratic  if  either 

a)  k  =  0     (vertical   intiidflnc©),   or 

b)  P^'^^*  P^^^  =  0  and  P^^V^"^  +  P^'^V^"^  =  0     . 

'  zx  xz  yx     zy  xy     yx  , 

When  this  happens,  we  have  h^   =  -  h^   and  h^   =  -h^   and  many  simpli- 
fications in  algehraic  tedium  result. 


-  8  - 


B-  Contimilty  Conditions  and  the  Coupling  Coefficients 

The  continuity  conditions  across  a  "boundary  require  that  the  tangentled 
components  of  "both  E  and  H  must  he  continuous  across  a  houndary  which,  in 
this  case,  is  taken  to  "be  the  plane  z  =  f  vrtilch  eeparates  the  l)   medium 


from  the  (u-l) 
yields 


th 


medium.  The  contiziuity  of  the  tangential  components  of 


where  we  have  suppressed  the  dependence  of  E  on  x  and  t.  It  is  understood 
that  (8)  is  to  hold  for  all  values  of  x  and  y  at  all  times.  We  are  now  in 
a  position  to  justify  our  assujnptlon  that  the  x-component  of  the  propagation 

vector  is  the  same  for  all  media.  E  is  a  function  of  x  only  through  the 

ikx  lk(w)x   ik'^~^/x 

exponential  factor  e   .  Hence  (8)  states  that  e     =  e        from  which 

we  see  that  k  '  =  k^~   =...=lc^   =k  since  k  is  defined  as  the  value 

of  the  x-conponent  of  the  propagation  vector  in  the  zeroth  medium. 

Because  we  are  dealing  only  with  the  electric  field,  the  two  continuity 

conditions  on  H  must  he  changed  to  additional  conditions  on  U  by  means  of 

Maxwell's  equations.  From 


dEz  . 

_  asv 

ay 

dz 

SEx. 

.  ^^z 

dz 

dx 

=   +iU)M.H 


=  +l«»)|j,H 


we  get,  realizing  that  the  electric  field  has  no  y  dependence  and  an  x  de- 
pendence equal  to  e 


<'>    »^  ifcn  -5 


as  (u-1) 


_  1 


t      ~J^     az 
S-i 


u-l 


h) 


U-l  ^    •- 


-iki 
C      « 
u-l 


"Ki] 
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Our  next  tatok  1b  to  substitute  (l)  into  eqaatlons  (8)  and  (9)  and 
thus  finally  express  the  continuity  conditions  on  the  electromagnetic 

field  as  a  set  of  linear  relations  whicH  give  the  amplitudes  C    In  terms 

(u-1)  ^ 

of  the  amplitudes  C^    ,  The  aforementioned  substitution  yields 

a)  y     c(-^>A("-^)e^^"'^«-l  =  t    C^'V^K'^^n^'h^i 
*— ,   n     nx  * — ^,  n   nx 

n=l  n=l 


^   '  '  ^— ,   n     ny  *— ,   n   ny 

n=l  n=l 


'   (i>-l)  '— ,  n    ny    n  [v)  *—.   n   ny   n 

[i.  n=l       *  R    n=l 

u     n=l  11    n=l 


where  we  have  defined  B   as 

n 


(11)  B^"^  =  h^"V"^-kA^"^ 
n     n   nx     nz 


It  should  be  noted  that  in  (lO)  the  exponential  e  '-"''    -■ ,  which 
appeairs  on  both  sides,  has  been  cemcelled.  The  solutions  of  equation 
(10 )  are,  in  matrix  form, 

(12)    5-  ^""^  «  ^M^^-l) 

where  C    and  C  '^~   are  tinderstood  to  be  the  four  element  column  vectors 
C^'^'^  and  C^*^-  .  The  components  of  f  ^'^^  are  given  by 

where 


(l^a)  ^ 


(») 
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\x    • 

^2x      • 

^3x     • 

^y    • 

^2y      • 

^3y    • 

4'^ 

•   ^2     ^2y 

•  ^3    ^3y 

^1    • 

^2        • 

33      , 

and  A    iB  ^        with  the  n   column  removed  and  replaced  "by  the  column 
zun 


»   _,r    t   (y-i)  m     my    '   (w-l)   m 


mx 


my 


St 


involving  the  properties  of  the  y-1   medium.  Thus,  for  example, 


(lU-b)  il 


2m 


\x 


H 


(u) 


A 


^(u-1) 
mx 

(u-1) 
my 


(u) 
3y 


^.^"^  .(u-l) 


B. 


(y) 


A 


»+y 


,(V-1),(U-1)    u(»>)a(u)     v(w)»(w) 


^"'^y''   "705=1)^  *''^;;~*'.  ^3'"^'3r  ^^    % 


B< 


(u) 
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C  -  The  Scattering  Matrix 

We  may  now  formally  carry  out  the  last  step  in  our  pro,=;rajn;  namely 
ohtaining  C^"  in  terms  of  the  C  .  The  connection,  in  matrix  form,  is 
simply  given  ty 

(15)  C<^^  =  1^"^  ^^'-"-^   ...  f^'^  =T^'^  C<°^  =  JJi  f^^^  C<°\ 

st 
Of  particular  interest  are  the  amplitudes  in  the  N  +  1   medium.   If  v;e 

suhstitute  N  +  1  for  U  in  (15),  we  will  get  the  amplitudes  of  the  four 

polarizations  in  the  homogeneous  isotropic  medium  ahove  the  anisotropic 

multilayer  in  terms  of  the  amplitudes  of  the  four  polarizations  in  the 

homogeneous  isotropic  mediiom  telow.  Because  of  the  special  nature  of  these 

two  media,  the  four  polarizations  wil  he  degenerate.  That  is,  equations 

(7)  will  hreak  down  and  yield  A^^  =  A   =  A   =m  •  0  so  that  the  ratios  of 

the  various  pairs  of  components  of  the  polarization  vectors  are  not  \inique- 

ly  determined,  as  is  well  known.  Hence  we  need  only  define  an:^'-  complete 

set  of  polarizations  in  these  media  such  as,  for  instance,  two  oppositely 

directed  circular  polarizations  or  two  perpendicular  linear  polarizations 

for  hoth  the  upgoing  and  downcoming  waves. 

Defining  f,  the  total  connection, or  coupling,  matrix  for  the  entire 
multilayer  by 

(16)  =T=|f  |C.).|(m)5(K)___=(2)=(l)  _   . 

where  care  must  he  taken  to  preserve  the  order  of  the  factors,  equation 
(15)  becomes,  finally, 

(17)  of*^>=i:  VO^"',    n=  1.2,3.'*. 

r=l 

Equation  (17)  determines  the  electric  field  above  the  multilayer  in 
terms  of  the  electric  field  below  the  multilayer  by  virtue  of  equation  (l) 
(if  v;e  define  any  complete  set  of  polarization  vectors  In  these  two 
isotropic  media).  Usiixilly,  however,  in  ionospheric  applications,  the 
amplitudes  of  the  two  downcoming  waves  in  the  zeroth  medium  and  of  the  two 
upgoins  waves  in  the  N  +  1^'''  medium  will  be  the  unknowns,  to  be  determined 
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in.  teiins  of  the  upgoing  waves  in  the  zeroth  medium  and  the  downcoming  waves 
in  the  H  +  1   medium.  The  amplitudes  of  the  do\/ncoming  waves  in  the  N  +  1 
medium  axe  usuclly  equal  to  zero,  for  in  most  cases,  there  is  no  wave  inci  - 
dent  on  the  ionosphere  from  atove.  At  any  rate,  it  is  desirahle  to  know 
the  amplitudes  of  the  outgoing  waves  in  terms  of  the  amplitudes  of  the  in- 
coming waves.  The  matrix  which  gives  this  connection  is  known  as  the 
scattering  matrix  and  is  generally  denoted  hy  S,  Mathematically,  S  is  de- 
fined hy 

■.(0) 
'3 

,(0) 

,(^^+l) 


(18) 


=   S 


,(I^1)J 


where  we  have  numhered  the  modes  so  that  the  subscripts  1  and  2  refer  to 
downcoming  waves  while  the  suhscripts  3  ^^^   ^  refer  to  upgoing  waves. 

It  would  "be  well  to  mention  that  in  the  case  of  the  multilayer,  the 
distinction  "between  the  ordinary  and  extraordinary  rays  in  each  layer  is 
purely  ar"bitrary.   It  is  only  in  the  continuously  varying  anisotropic  medium 
at  relatively  high  frequencies  that  the  concept  of  a  ray  is  at  all  sensihle; 
hence  the  notation  of  ordinary  and  extraordinary  rays  cannot  "be  applied  to 
the  multilayer  to  advantage. 

It  remains  now  only  to  give  the  components  of  the  scattering  matrix 
S  in  terms  of  the  connection,  or  coupling,  matrix  T.  This  can  be  done  most 
simply  "by  solving  the  four  linear  simultaneous  equations  given  by  (17) 
fo;c(?)4^'\  Cf  ^).  Ci^*^^   in  ter.s  of  C^^J  C^^^.  o[^^'\    cf*^^ 
NOW,  C ^^^  '  =  f  C^  ,  or  in  partially  expanded  form. 


(19r:) 


+  t 


where  t  ,  t  ,  t  ,  t  are  the  two-hy-two  matrices 
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(191)) 


'V^21  •   ^22/ 


Solving  (19)  for  the  outgoing  wares  in  tenns  of  the  incoming  ones  yields 
the  following  expression  for  S 


(20) 


,(0) 

1 

(0) 

2 

(H+1) 


=  § 


a    D 

d  c  a   D 


if   t      is   non-singuler. 

Taking  the  prnplitndes  C^  '      ^  and  Cp  of  the  downcoming  waves   in  the  N+l 

medi-om  to  he   ^ero   (which  U8i:ially  is   the  case),   equation   (20)  simplifies  to 


(21a) 


.(0)' 
,(0) 


a       h 


which  gives  the  amplitudes  c|  and  C^   of  the  reflected  waves,  and  the  amplitudes 

C^^'^^and  CP"*"   of  the  transmitted  waves  in  terms  of  the  amplitudes  c).     and  Cj^ 

of  the  incident  waves.   Since  the  media  above  and  helow  the  multilayer  are  isotropic, 


we  may  define 


.  (N+l)_  .  (W-l)_  .(N+1) 
^d    -  ^1    -  h 


u(0)_  ,,(0)_  .  (0) 
h_  =  h^  -  h^ 


u 


,  (N+l)_  A^+l).   v,(N+l) 
h.    -  h^    -  h^ 


u 


and 


E 


in 


=  (c'°)A(<"«i°)li<"). 


i(-ujt+kx+h^°^z) 
^       u 

(0). 


for  z  <  0 


(21h) 


for  z  <  0 


,^\    /^\  /r.s    /^\   i (-"»t+icx+h]  ^z) 


So  that  eanPvtlons  (21a)  yield,  as  our  final  result,  the  connection  "between  the  inci- 
dent field   E   end  the  reflected  and  transmitted  fields  E^^^  and  E^^, 
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D  -  Application  to  the  Ionosphere 

In  section  C,  the  coupling  and  scattering  matrices  were  derived 
for  an  arbitrary  anisotropic  multilayer  subject  only  to  the  restriction 
that  ^1  was  a   scalar.   In  practice,  one  of  the  most  important  anisotropic 
media  considered  is  the  ionosphere,  in  which  the  conductivity  tensor  W' 
is  due  to  free  electrons  in  a  damping  medium  permeated  hy  a  fixed  uni- 
form magnetic  field  and  in  which  p,  and  e  are  scalar  constants.  F  is 
found  hy  solving,  first,  for  the  velocity  v  from  the  equation  of  motion  (in 
M.E.S.  units) 

i     -      —  -iu)t      —    — 
mv+  m»v=  -eEe     -iJievx  H 


of  an  electron  in  a  resistive  medium  under  the  influence  of  a  fixed  magnetic 

-iujt 


field  H  and  an  oscillating  electric  field  E  e 


Here  m  is  the  mass  of 
the  electron,  v   is  the  collision  frequency  and  e  is  the  magnitude  of  the 
electronic  charge  (the  sign  having  already  been  taken  into  account). 


Using  the  relation 

=  -  -i<ut     —      — 
CTEe     =  j  =  -Ne  V 

0 

—  ,  —  -ittit 

for  the  current  density  j  in  terms  of  the  impressed  electric  field  E  e 

where  N  la  the  number  of  electrons  per  unit  volume,  f- can  be  obtained  from  v. 

Letting  w,  UJ-  and  W-,  represent  the  angular  velocities  associated  respect- 
H      Cr 

Ively  with  the  frequency  of  the  incident  wave,  the  gyromagnetic  frequency 

and  the  critical  frequency,  i>5. 


^  ue  IHI 
m 


}      Cr 


~    J     em         ' 


(22)  uj„ 

^   '  H       m    y   Cr    V   cm 

and  letting  a,  p,  V  be  the  direction  cosines  of  the  magnetic  field  H  ,  r 

is  given  by 

2 


ic« 


0-  = 


Cr 


((^iu)[(^+iu)^-«"„^J 


(«M-iu)2-ui/a^ 


•   -idu+iyiugY   -uig   00 


i(uH-iy)uj   V_Ui  ^Pa;    ("J+iu)^-Ui--^^ 


H 


E 


H 


-i(tt)+iu)o)gP-ujg  Va  ;   i  (uj+iu  )ujga -Uig     -yp; 


(23) 


-i(uj+iy)ojgar-"^^p7 
((U+iy)^-  «^^V^ 


U 
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With  <r  given,  p  is  determined  "by  equation  (^b): 

{2k)  P  =  ".^|x«  (!  +  ^  (T  ) 

for  scalar  c.  From  P,  the  scattering  matrix  for  the  ionospheric  multilayer 
can  he  computed. 

That  the  multilayer  can  he  used  to  approximate  a  continuously  varying 
medium  follows  from  the  fact  the  scattering  matrix  is  a  continuous  functional 
of  the  cliarge  density  and  collision  frequency.   It  is  a  consequence  of  this 
fact  that  as  the  number  of  layers  in  the  multilayer  tends  toward  infinity  so 
that  the  multilayer  approRChea  a  continuous  mediijm,  the  scattering  matrix 
for  the  multilayer  approaches  the  scattering  matrix  for  that  continuous  medium. 
Hence  we  are  led  to  expect  that  for  a  reasonably  Inrge  number  of  thin  layers, 
the  scattering  matrix  for  the  multilayer  will  he  a  good  approximation  to  the 
scattering  matrix  for  the  continuous  medium  it  represents.  This  will  he  the 
case  if  ench  layer  is  hut  a  small  fraction  of  a  wavelength  in  thickness. 
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3.  THE  CONTIKUOUSLY  VARYING  MEDIUM  AT  KIGE  FKEXjUEITCIES 

The  multilayer  provides  em  approximation  to  the  actual  scattering 
matrix  for  the  continuous  medium  in  the  low  frequency  range  which  cannot, 
in  general,  he  treated  hy  other  methods.   In  the  hi^  frequency  range, 
however,  a  solution  in  series  form  exists  which  is  essentially  of  the 
W.K.3.  type.   The  n   term  in  the  series  is  of  the  order  of    «    as 
«u  — >oo,  as  will  he  seen  subsequently.   As  a  result,  it  is  possihle  to 
effect  a  considerahle  simplification  in  the  analysis  hy  expanding  the 
quantities  appearing  in  each  term  in  the  series  in  inverse  powers  of  '^  » 
The  importance  of  such  a  simplification  can  he  appreciated  when  It  is 
realized  that  the  quantities  appearing  in  the  original  series  involve 
the  solutions  of  a  quartic  equation  and  vtrious  complicated  functions 
of  these  solutions.  On  the  other  hand,  hy  virtue  of  the  expansion  in 
inverse  powers  of  ui      it  is  possible  to  ohtain  relatively  simple  useable 
series  expansions  for  the  total  coupling  and  scattering  matrices  as  ex- 
plicit functions  of  all  ionospheric  parameters. 

Maxwell's  equations  can  he  transformed  into  a  system  of  four  coupled, 
linear  first  order  differential  equations*  in  four  amplitude  functions 

^i,(z)  by  substituting 

k  _  i[-«<t+kx+/^h   (C)dC] 

,       E(x,2.t)  =  Yl    V   (z)iA   (7)e  °    ^ 

(25)  n=l    ^  ^ 

k  _  i[-U)t+kx+/\   Odt] 

E(x,7,t)  =  Y.    V^(z)D„(z)e 
n=l     "         ^ 

into  Maxwell's  equations,  v;here  A  ,   D     and  h     satisfy  the  algebraic 
^  n'      n  n 

relations 

(26)  (^)  -Vny  =  -"^\x  ^'^  -Wv  =  k  ^^'^r?.__ 

(c)  kA^y  =  -"i.D^^  (f)  kD^  =  ^  (PA^), 

which  are  equivalent  to  Maxwell's  equations  for  a  plane  wave 

i  [-u)t+k^+h^z] 

E     =  A       e 
n         n 


^^^^  if-cut+k  +h  z] 

_         —  *-  X     n  -^ 


H     =    D 
n  n 


*H.  Bremmer,    Loc.   cit. 
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traveling  in  a  homogeneous  medium  whose  properties  are  the  same  as  those  of 
the  actual  medi\m  at  the  point  z.  This  suhstitution  yields,  for  n  =  l,2,3.^t 
the  differential  equations  satisfied  by   ^  (z) 


dV 


(28)      !^ 


dz 


=  -e 


-•"  "     r  a  (z)V^(z)e  °  "^ 


■'o  n 


m=l 


where 
(29a) 

(29h) 


A 

mn 

mn        A 

^l,' 

^2x* 

%• 

^. 

A   = 

^17' 

^ay- 

*3.- 

^y 

hhr 

V2y 

"3*3^ 

Vi.y 

h' 

»2- 

^3- 

\ 

and  A   is  equal  to  A  with  the  n   column  removed  and  replaced  "by  the 


column   dA 


mx 


dA 


my 


dB 


d7  •  dz  •  dz^  m  my''*  dz 


with 


(29c) 


B  =  h  A   -  kA 


"n    n  m 
Thus,  for  example 


(29d) 


a 


"23  "  A 


nz 


h' 


^3x 
dz   • 

d- 


A 


33C' 


^37- 


fc(V3y^'  V3i 


dB- 
dz 


B«,  t 


^/^x 


V/.y 


B, 


The  a   given  in  (29)  are  exactly  the  seme  as  those  given  hy  Bremmer.  However, 

nm 
they  are  here  written  only  in  terms  of  the  electric  field,  the  tangential  com- 
ponents of  the  magnetic  field  D^  and  D   having  heen  eliminated  with  the 
help  of  (26a  and  h). 

The  solution  to  (28)  is 


(30)   \i^)   =  E  Mnm^^)  \^Q) 
m=l 
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where  M  (?)  is  given  by  the  series 

00 

(31)     M„^(z)  =  ri     "A (2)       .  where 
nin  r~"       nn  ' 

J=o 


Tifii  zun 


nin 


nm  s    =1       9       =1  ' 


(32) 


J-1 


•expi^^^h^Od:  *  |:n^^^  \^^)^^  *Jo\^^^^^l    '°"  "5  > 

That  (?0)  is  formally  a  solution  of  (28)  cen  "be  verified  hy  airect  substitution. 

The  series  (31)  for  M  (z)  is  given  in  terms  of  the  h  and  the  a  „. 
*   '      nm »  '    ■=  r         rs 

The  former  are  the  solutioas  of  a  quartic  eq\jation;  the  components  of  the 

matrix  a  are  given  in  (29)  as  the  ratio  of  two  four  by  four  determinants 

whose  elements,  in  turn,  are  functions  of  the  h  .   It  is  thus  evident  that 

the  formal  solution  (30)  is  rather  impractical  for  demonstrating  the  overall 

dependence  of  the  scattering  mstrix  on  the  ionospheric  properties. 

On  the  other  hand,  if  p  is  expanded  in  terms  of  inverse  powers  of  the 

frequency,  h  and  a  can  also  be  obtained  as  series  expansions.   It  is  then 
n 

possible  to  obtain  relatively  simple,  useable  series  expansions  for  the 
total  coupling  and  scattering  matrices  as  explicit  functions  of  all  iono- 
spheric paramters.  These  expansions  will  converge  at  freq.uencies  greater 
than  both  the  gyromagnetic  and  critical  frequencies,  which  is  exactly  the 
frequency  range  in  which  (31)  is  rapidly  convergent. 

Thus,  writing  ^  as  a  power  series  in  the  dimensioniess  parameter 
\  =  <*w/"',  we  have 

(^  A.        H  r=o 
i  being  the  unit  matrix.  The  first  few  terms  of  this  expansion  are' 


(3H) 


P  = 
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2                tt,     2 

+a     ;         +V 

-V    ;         •»» 

+  P    ;         -a. 

;     -3 
•      +a 

u.    2 

a^^a^-l     : 
ap-2aV     ; 

aV+2aP       ; 

Pat2aY  ; 
a^.p^.i     . 

pV-2aa       ; 

Var 2a p 
Vp+2aa 

a^-HV^-l 

a(3-3a^-a^)    :  V-3aV3apa  ;  -p+3a  p-3ava 


2  2. 


-V+3a^-3aap   ;  a(3-3P  -a  )  ;  a-3a  a  -3aVP 


2  2. 


P«3a^P  -3aar('  ;  -<i  +3a  a-3apV  ;  a(3-3v  -a  ) 


+  •♦•  + 


where  a  is  defined  by  a  =  ;^.  Substituting  this  expression  into  the 

H  ^ 

three  equations  (5),  we  can  solve  for  h^  and  A^  as  series  expansions  in  \: 


(35) 


^  ^)  ^  =  -\^°°«^*^„^n 


^       (r).(r)_5  ^leose.E^^^^^^ 
-'n  OX  * —  " 


00 


r=0 


r=0 


1))  X  =  Fa 


7(rKr 


r=0 


U) 


^=5=1  and  f  ^  =  5  K  =  -1  80  thst  the  term  -  ^  -r  cos  Q  in  (35) 
expresses  the  fact  that  as  o>— >  oo  and  hence  X  — ^>  0,  h^ — ^>  -  ^^-^  cos  9, 
the  values  of  the  z-component  of  the  propagation  vector  in  free  space  for 
the  downcoming  and  upgoing  waves  respectively. 

Carrying  out  the  indicated  suhstitutions  and  suhsequently  equating 
the  respective  coefficients  of  X  to  zero  in  each  of  the  three  equations 
(5)  yields,  for  each  r,  a  set  of  three  simultaneous  equations  which  can  he 

solved  for  h^'"'*'  \  a  relation  Vetween  A^    'and   A^f  ' ,   and  finally, 

n  nx       ny   ^(r_i) 

for  r  >  1,  a  relation  between  the  three  components  of  A     .  Thus, 

after  p.  little  collecting,  we  find  that  there  nre  two  equations  con- 
necting; the  three  comr.onents  of  A^^',  while  hj;   is  uniquely  deter- 
«=  n  n 

mined.   This  corresponds  to  the  fact  that  equations  (5)  det^-'rinine  the 

h  unia-oely.  while  vieldirw?  only  two  relations  between  the  three  com- 

n  ■  "    _ 
ponents  of  A  .  (i.e.  the  normalir.etion  of  the  polarization  vectors 

A  remains  arbitrary.)  Although,  at  first  sight,  it  might  seem  strange 
that  the  two  equations  connecting  the  components  of  A^*^  are  obtained  from 
different  orders  in  the  expansion  of  equations  (5)f  a  little  consideration 
makes  it  clepr  the t  this  lack  of  symmetry  is  necessitated  by  the  require- 
ment that  if  the  frequency  of  the  incident  wave  becomes  infinite  or  if  the 
charge  density  vanishes,  the  polarization  vectors  A  must  become  perpen- 
dicular to  the  propagation  vectors.  That  is,  lcA^y+  \^nz~  0  ^^  ^  =  0, 


(r)      (r) 
or  if  N  -   0,   Thus,  the  equrtion  connecting  A    and  A    stems  from  a 

different  source  then  the  other  equation  connecting  the  three  components 

of  A   ,  and  we  must  not  "be  surprised  to  find  this  lack  of  synimetry. 


The  exoansions  for  h  and  A  have  been  evaluated  to  the  order  \ 

n      n 


h  =-^? 


«)_ 


H 


\  ^n  C 


cos 


G  +  \| 


Cr 


n  ZttJ-CcoeO 

n 


-y}x 


u)„  2 

Cr 


n  2U'_Cco 


5 


A   =  ?  cos©  -  X  ? 
nx    n  n 


u)  2 

Cr 


36) 


2uj„  cos  6 


ny 


_  ..    X  i -S ^  ,  2  iCn  K  '^     4  it  t^    n  Cr  f 

^n  ~   2  cos  0       2  I  itco620     ^n  U)  co80  "  oTT"  f 

^n         V  '^n       H    '^n  H   ^ 

2 


A   =  sin  e  +  OX  +  OX 
nz 

where  we  have  so  normalized  A  that  A^  '=  A   =  0.  The  notation  used  above  is 

n      nz    nz 

summarized  on  Table  I. 


S,=  l    §2=^   5  3=-l 
Xl=l    ^2  =  -^   ^3" 

:,  =  !       :3  =  -i     C3—1 


! 


=  -1 


Ul 


X  = 


H 


Ul 


.  «)„  (z)  = 
'  Cr     > 


i^(z) 


^   i?(z)  is  the  collision  frequency. 


'"  '    h.a(»' 


a  =  a^  cosO  +  Pi  C  +  "VsinO  = 
n    ''  n         n. 


H 


cos  i     -   asin«  -  "V  ^  cose 
TABLE  I 


It  will  be  noticed  that  0  is  the  angle  made  by  the  incident  wave  normal 

of  the  n   mode  with  the  direction  of  the  magnetic  field. 

The  expansions  for  h  A   and  B  are  simply  obtained  from  those  given 

n  ny  2    ^ 
"by  (36).  They  are,  to  order  A  : 


u)„  U)_fl.  cos 

nV  =  -  X  ^^  -  °°^  ^  -  ^  >n  20  cos  i 


cos  9 


-^^^ 


(37) 


uj,,cos  e  ,  - 

K     \  _ii +  i^  iL 

f    2j      n  <D, 


2C 


W, 


B  =  h  A  -k.^   =   1  !^  +  X  -^# 
n    n  lut   nz     X  C      w  c 

H 


cos 


U)  2      2  •^ 
cr  1+cos  t    / 

(11.^2  C032  9  I 


Having,  now.  the  exoansions  for  A  ,  A  ,  h  A   and  B  ,  we  csin  compute 
°»    »      '■  nx   ny   n  ny      n 

the  expansions  for  ^nm.  Zik  .  and  hence  a   .  Orfing  to  the  fact  thst  the 
^  '  lun 

first  two  terms  in  the  expansions  for  A^,  A  ,  h^A   and  B^  are  independent 
of  z,  the  series  representing  a   will  begin  v/ith  the  term  in  X^  ^^  \,\xbX 
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00    /  \ 


(T8)   a  (z)  =  5~  a^^^(z)\' 


Denoting  the  err:ionential  phase  fnctor  for  the  Jth  term  in  (32)  by 

(39)     ^(o)  ^  1 
nm 

(J)  it-io\(OdC  +  ^  h^  (c)dc  *  ^X^^)'^^! 


$ 


(^0) 


n.9^» . . .  fS  1   2^»ni 


=  e  k=l       k 


the  Buhscripts  n,  s^^,  •••»^j«i»  "^  taklpg  on  the  values  1,2,3,^,  we  have 
for  }^H^) 


nm 


nm      nm 


nm^'     Jo   l*nm^l'  *— ^  nm  ^  1 ' 
'  r=2 

M^J^z)  =  (-1)''  7"  ...7"    \   dz.  \  dz,  ...  )  •"  dz,$  , 

nm  '  '        g^r^    8   =1  '°  Jlnj^8^,...,8  ^^,m(z,Zj^,... 

r«2  '^^l  ^      r=2  ^1^2  "^        r=2  ^j-f  •' 

In  (40)  we  have  not  expanded  the  exponential  phase  factor  in  powers  of  X, 
■because  the  coefficients  in  such  an  expansion,  involving  integrals  of  the 
propagation  vector,  will  he  roughly  proportional  to  the  thickness  of  the 
ionosphere  in  wavelengths.  At  high  frequencies,  this  will  be  very  larg:e, 

so  that  the  coefficient  of  X^  in  the  expansions  of  the  ^  's  will  he 

r  — 

very  much  larger  than  the  coefficient  of  X   in  the  expansion  of  a  .  Leaving 

the  exponential  unexpanded  will,  on  the  other  hand,  insure  that  terms  of 

different  order  in  X,  in  the  expansion  of  M  (z),  will  he  of  different  order 

'  nm   ' 

of  magnitude  in  cases  of  practicad  interest. 

It  can  "be  seen  from  (UO)  that  M  ^     involves  the  product  of  j  a's  • 

nm  "  2 

thus,  since  the  series  expansion  for  a  begins  with  the  terms  in  X  ,  the 

expansion  for  M  "^  begins  with  the  term  in  X   as  mentioned  earlier.   Con- 
sequently, the  contribution  due  to  multiple  reflections  decreases  with 

increasing  frequency, 

(2) 
The  leading  term  a   (z)  in  the  expansion  of  a   has  been  calculated: 
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2  v»  .»    ,  ..  .2 


Cki  t  (2)/    V  Cr  (^       ^nlm     n  m  SnSm     ^nXm   ^         ,  n 

(Ul)       a^  (z)  =  -  j— ^—     (1  +  2  +       ,  ^     2,     ^  *>"  SSTT^ 


H  2  cos  e  H       '^m 


where  oi         and  v  are  functions  of  z  and  the  prime  denotes   differentiation, 

2 
From  this  we  ottain  M     (z)  to  order  X     "by  adding  the  contributions   due  to 
/  /\  \  /  ^  \  Dm 

M       (z)  and  M       which  ere  the  only  terms  that  contribute  to  order  X 

,,..   M     (z)  =  6     -X^pdze'^'     ^  ■"  Jf^(zJ      . 

v'+c^     nmnm/oi  nml 

In  (U2)  the  factor  e~^/o^^^^  ^^'^^  was  taken  into  the  integral 

.i/X(Odt    i[r\(0d5*/^\(0d:]      i/^i[-h^(o*h^(:)Jdt. 

e  •  e    -1-  =  e 

The  scattering  matrix  for  a  continously  varying  anisotropic  medi-um  of 

2 
thickness  L  is  now  ohtainahle  (to  order  X    )  from  M  (L).  -l^efining  the  amplitudes 

C^(z)  hy   the  relation       ^^   ^^^^        ^^^   (^j^j 

{U3)  C^(z)e     ^  =  \(^)«         " 

we  have,    from   (25),  _ 

i4  _  i[-u)t+kx+h   (z)zj 

(i:i')  E(x,z.t)  =  51  C    (z)A    (z)e  "" 


n=l 


The  coupling  matrix  T  for  the  medium  is  defined  "by 
h 

(^5)     C  (I.)  =  Yl     T„„C  (0)  . 
n      ^~,   nm  m 

m=l 

Substitution  of  (U3)  into  (30)  and  setting  z  =  L  yields 

-ih  (L)L  iA  (Od? 

(i*6)  T   =  e   ^    e  °  "     M  (L) 
nm  nm 

The  scattering  matrix  S  is  given  in  terms  of  the  elements  of  f  hy   equations  (19h) 

and  (20).   If  the  incoming  field  is  Incident  from  helow,  the  reflected  and  trans- 

-  2 

mitted  fields  sre  given  "by  equPtions  (21),  The  components  of  f  to  order  \  are 

reafMly  o^htained  hy  substituting  (^2),  (^1)  and  the  first  of  equations  (36)  into 
(i|6). 
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No  special  caees  will  "be   considered  here  "but  it  can  "be  concluded  in 
general  from  the  a"bove  results  that  there  is  an  anomalous  effect  when  the 

direction  of  either  an  incident  or  reflected  wave  coincides  with  the  direction 
of  the  earth's  magnetic  field.  Attemiation  in  the  medium  is  roughly  proportion- 
al to  the  path  length  of  a  geometrical  optics  ray  in  the  medium  and  depenos  on 

the  product  of  the  collision  frequency  and  charge  density.  Strictly  speaking, 

2 
the  absorption  factor,  to  order  \   ,  is  Just  the  ahsorption  that  would  he  suffered 

"by  a  geometrical  optica  ray  which  is  incident  at  angle  6  and  traverses  the  medium 

tmdevlated  except  for  a  transition,  at  the  height  z, ,  from  the  m   mode  to  the 

n  ,  provided  that  the  m   mode  is  upgoing  and  the  n   is  downcoming.  Higher 

order  terms  take  into  consideration  the  curve ture  of  the  path  in  the  medium  and 

the  effect  of  the  earth's  magnetic  field. 

ProB  (^1),  it  is  seen  that  even  aside  from  the  attenuation  effects  Just 

mentioned,  the  scattering  matrix  will  depend  on  the  collision  frequency  through 

the  term  .   2 

tl2  ,    n 

^         8«i)„co8«5 
H   '^m 

and  that  this  effect  will  depend  also  upon  the  magnitude  of  the  earth's  magnetic 

field  and  the  direction  of  the  incident  wave  with  respect  to  this  field.   Carrying 

out  the  expansions  to  higher  order  will,  of  course,  yield  better  acciirecy.  However, 

2 
the  major  effects  already  appear  la  the  \   term  in  the  frequency  expansion. 
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